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ABSTRACT 

Errors in the kinematic distances, under the assumption of circular gas orbits, were estimated by 
performing synthetic observations of a model disk galaxy. It was found that the error is < 0.5 kpc 
for most of the disk when the measured rotation curve was used, but larger if the real rotation curve 
is applied. In both cases, the error is significantly larger at the positions of the spiral arms. The 
error structure is such that, when kinematic distances are used to develope a picture of the large scale 
density distribution, the most significant features of the numerical model are significantly distorted 
or absent, while spurious structure appears. By considering the full velocity field in the calculation of 
the kinematic distances, most of the original density structures can be recovered. 
Subject headings: Galaxy: disk — Galaxy: kinematics and dynamics — Galaxy: structure — ISM: 
kinematics and dynamics — MHD 



1. INTRODUCTION 

Since the classic work bv lOort et al.l l)1958|) . there have 
been many attempts to use the kinematic properties of 
the diffuse gas to determine the large scale spiral struc- 
ture of the Milky Way. Very early in the study of the 
Galaxy, it was determined that the orbits of the disk 
components of the Galaxy were not very different from 
circular, with an orbital frequency that decreases mono- 
tonically as a function of galactocentric radius. These 
facts allow the use of the kinematic distance method as a 
first approximation to map the gaseous component of the 
galactic disk. Two of the main strengths of this method 
(that it can be used for a very large fraction of the Galaxy 
and that it can be applied to the gaseous component of 
disk, which is notoriously difficult to obtain a distance 
to) make it particularly useful for this goal. Neverthe- 
less, it was soon realized that the deviations from circu- 
lar orbits, however small in absolute value, might have a 
strong impact on how we see the Galaxy. 

One of the first difficulties of the kinematic distance 
method appeared in the determination of the rotation 
curve, namely, in the fact that the circular rotation ve- 
locity measured for positive galactic longitudes (north- 
ern Galaxy) did not match the one measured for nega- 
tive longitudes (southern Galaxy). The simplest way to 
reconcile these observed rotation laws is to take their av- 
erage, assuming that the differences gener ated by non- 
axisymmetric structure will cancel out. iKerrl l|1962f) 
showed that this approximation leads to large north- 
south asymmetries that, given their heliocentric nature, 
seemed unlikely. It became clear that the complex kine- 
matic structure revealed in the diffuse gas surveys, like 
the presence of gas at forbidden velocities, or the oscil- 
lations in the rotation curve, was itself a consequence of 
the spiral structure that was being sought. Given the im- 
portance of the rotation curve for the understanding of 
the galactic dynamics, in addition of the determination 
of kinematic distances, a different approach to measuring 
the rotation curve was needed. 

A frequently used method to obtain the galactic ro- 
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tation curve involves measuring the full velocity field 
of discreet sources that might share the motion of the 
diffuse gas (young objects like H II regions, for exam- 
ple), and then averaging the so obtained azimuthal ve- 
locities. Thi s approach has lead to models of the ro- 
tation curve ijBrand fc Blitzlll993t iMaciel fe Lagdl2005|) 
that might more closely trace the real mass distribu- 
tion of the Galaxy, but introduce new sources of error 
when used to obtain kinematic distances. Neverthe- 
less, generally accepted mode l s of the spiral structure 
l|Georgelin fc Georgelinl Il976t iTavlor fc Cordesl Il99^ 
have been obtained using th is assumpt ion. (For a nice 
review of the early work, see lKern fl969) 

Another approach involves the modelling of the non- 
circular motions of the gas, instead of forcing the 
assumption of circular orb its. In a recent paper, 
iFoster fc MacWilliamsl l)2006j) used an analytic approach 
for the velocity field of the outer Galaxy. In this work, 
a numerical model of the galactic disk, with full MHD, 
is used to further explore the effects of non-circular mo- 
tions in the image one would obtain of the Galaxy when 
we rely on the kinematic method for distances. An ob- 
server is imagined inside the numerical model, which is 
assumed similar to the Milky Way, and the analysis that 
this observer would perform is reproduced. In a brief 
description of the numerical simulation is presented; in 
the selection of the observer's position is described, 
and how the measurement of the rotation curve was em- 
ulated; 21 presents an analysis of the errors in the kine- 
matic distances and how they affect the image the ob- 
server generates of his/her home galaxy; finally, [^sum- 
marizes the results. 

2. THE SIMULATION 

The numerical model used here is d escribed else- 
where ijMartos et al.ll2004alri lYanezll2005|) . and only an 
overview is presented here. 

The initial setup consisted of a gaseous disk with an 
exponential density profile in the radial direction, with 
a scale length of 4 kpc. The density at the position of 
the Sun (at 8 kpc from the galactic center) was taken 
to be 1.11 cm~ 3 . The equation of state for the gas was 
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isothermal with T = 10 4o K. This disk was threaded by 
an azimuthal magnetic field, with strength given by the 
relation, 



Pb =Pm- 



(1) 



where ps is the magnetic pressure, n is the gas density, 
Pm = 1-43 x 10~ 12 dyncm~ 2 , and n c — 0.04cm -3 . Equa- 
tion (JTJ yields B — 5.89 /iG for r = Rq. This intensity, 
and the magnetic field geometry, were set only as initial 
conditions, and were allowed to evolve in time according 
to the ideal MHD equations. 

The gas initially follows circular orbits, with a velocity 
given by the equilibrium between the background grav- 
itational potential, the thermal and magnetic pressures, 
and the magnetic tension: 
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where is the azimuthal velocity, p = m e ffn is the gas 
mass density, m e g = 1.27rn# is the mean particle mass, 
Pt is the thermal pre ssure, and $ is the grav itational 
potential described bv iDehnen & Binnevl (jl998t model 

The equilibrium was then perturbed by the two- armed 
spiral potential described in IPichardo et alJ ((2003). The 



simulation was performed in the perturbation reference 
frame, which rotates with an angular speed Qp = 
20kms~ 1 kpc~ 1 l|Martos et all l2004albT ). It is worth 
mentioning that the perturbing potential was calculated 
as a superposition of oblate spheroidals, and so it does 
not have the usual sinusoidal profile. Also, the param- 
eters that describe the perturbation (total mass in the 
arms, pitch angle, pattern speed, etc.) were constrained 
by those authors so that the pattern is self-consistent in 
the stellar orbits sense. 

The MHD equations were solv ed using a version of 
the ZEUS code (|Stone fe Normadll992a1bh . a finite dif- 
ference, time explicit, operator split, eulerian code for 
ideal magnetohydrodynamics. The numerical domain 
consisted in a two-dimensional grid in cylindrical geom- 
etry, with 500 2 points. The numerical domain extended 
from 1 kpc through 15 kpc in radius, and spanned a full 
circle in azimuthal angle. The boundary conditions were 
reflecting in the radial direction. 

Figure shows the simulation after 1 Gyr of evolution. 
The most important characteristics of the simulation at 
this stage are: 1) although the perturbation consists of 
two spiral arms, the gas forms four arms (two pairs with 
pitch angles of 9° and 13° each, as opposed to the per- 
turbation with a pitch angle of 15.5°); 2) a high density 
ring is formed at r ks 4 kpc; and, 3) a low density ring 
is formed near corotation, at r « 11 kpc. Again, details 
of the simulation and the physical phenomena related to 
these structures are discussed elsewhere. The corotation 
low density ring was found (using a different background 
potential) by Martos & Yanez (private communication). 
A study of the neccesary conditions for the formation of 
such ring, its physics and consequences for our Galaxy is 
p resented in Mar t os (2006 , in preparation). 

iGomez fc Coxl l)2004albf) also performed large scale 
simulations of the Galaxy. Since their numerical model 




Fig. 1. — Density distribution of the simulation after 1 Gyr of 
evolution. The grayscale shows the gas density in units of cm -3 . 
It is noticeable that, although the perturbation has only two arms, 
the gas response shows four arms. The Sun symbol shows the 
position chosen for the imaginary observer. 



was three-dimensional, they were able to study some phe- 
nomena (like the vertical motions associated to the hy- 
draulic jump behavior of the gas near the spiral arms) 
that could have an impact in the dynamics of the gas 
near the midplane. Nevertheless, the focus of their model 
was to study those phenomena, and an emulation of the 
Milky Way Galaxy was not a priority. Specifically, their 
three-dimensional numerical grid, a necessity in their 
work, restricted the spatial resolution achivable in the 
midplane. This, together with the low value for fip used, 
did not allow the formation of four spiral arms as a re- 
sponse to a two-arm spiral potential (in order to obtain 
four gaseous arms, their model included a four-arm po- 
tential). In the present work, the model was restricted 
to the galactic plane so that sufficient resolution can be 
reached. 

3. THE SYNTHETIC OBSERVATIONS 

Local maxima in column density (N) vs. galactic lon- 
gitude (I) plots for the diffuse gas are usually interpreted 
as the directions at which the line-of-sight is tangent to 
spiral arms. The N vs. I distribution an imaginary ob- 
server would see if placed inside this model galaxy is 
presented in Figure [21 By moving the observer around 
the solar circle, at 8 kpc in the numerical model, the 
number and positions of the local maxima can be fit- 
ted to the observed values for the diffuse gas. In this 
work, the chosen directions were thos e tangent to the 
locus of the spiral arms proposed by iTavlor fc Gordesl 
(1993), namely, I = -76.5°, -51.6°, -32.4°, 32.7° and 
47.1°. It was found that by choosing the position shown 
in Figure ^ f° r the imaginary observer all but one of the 
column density local maxima in the model fall within 
3° of these qu oted directions. If the t angent direc- 
tions quoted by Drimmcl & Spcrgcl (2001y are adopted, 
namely I w -80°, ±50° and ±30°, all but one of the tan- 
gent directions yield an even better fit. (The difference 
between the ill-fitting tangent in the model, at I w —72°, 
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Fig. 2. — Column density, in arbitrary units, vs. galactic longi- 
tude in the simulation. By moving the imaginary observer around 
a circle of radius Rq = 8 kpc, the local maxima can be matched to 
the observed directions tangent to the spiral arms (thick vertical 
lines). The dashed area at |/| < 7° shows the simulation inner 
radial boundary. 
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Fig. 3. — Line-of-sight component of the velocity field as a 
function of galactic longitude and (real) distance to the observer, 
with contours every 50 km s — 1 (the thicker contour marks the u r = 
0km s _1 level). The dotted line shows the locus of the tangent 
points, while the circles show the positions at which the terminal 
velocity is reached. 



and the quoted direction is in fact smaller than the width 
of the feature obs erved in 240 /xm. See, fo r example, 
lDrimmelll2000j and lDrimmel fc SpergefeOOll l 

3.1. The rotation curve 

Once a position for the observer is chosen, the next step 
toward calculating the kinematic distances is to adopt 
a rotation curve for the simulated galaxy. For the in- 
ner galaxy (r < Rq), the standard procedure consists in 
searching for the terminal velocity of the gas, i.e., the 
maximum line-of-sight component of the velocity (mini- 
mum, for negative longitudes). If one assumes that the 
gas orbits are circular, the terminal velocity arises from 
the point at which the line-of-sight is tangent to the or- 
bit, and so, the galactocentric radius of the emitting gas 
is known. Under this assumption, the circular rotation 
curve for the galaxy is given by: 

v c (r) = v t (l) + v Q sin(7), (3) 

where v c is the circular velocity, v t (l) is the terminal 
velocity for a given galactic longitude, v Q is the velocity 
at the solar circle, and r = Rq sm(l) is the galactocentric 
radius of the tangent point. 

At this point, a choice between two options for the 
value of the circular velocity at the solar circle has to 
be made. One option is to take the circular velocity 
consistent with the background gravitational potential 
(vq = 220 km s _1 ). This option has the disadvantage 
that the gas in the evolved simulation will stream by the 
observer [although this is not necessarily wrong, since 
the presence of gas at f orbidden veloc i ties in the I — v 
diagram is well known l|Linbladl Il967t iBlitz fc Soergell 
1991)]. Nevertheless, it was decided in this work to 
take a second option, which is to take the value for vq 
(= 225km s" 1 ) given by the azimuthal velocity of the gas 
in the evolved simulation at the position assigned to the 



imaginary observer, since this choice would more closely 
mimic the procedure used to deter mine the Local Stan- 
dard of Rest from galactic sources l|Binnev fe Merrifieldl 
1998). There will still be streaming gas, but this will hap- 
pen in the radial directio n only. Such radially s treaming 
gas has been reported bv lBrand fc Blitzl l)1993[K 

Figure |21 shows the line-of-sight component of the ve- 
locity field. The Figure also shows the actual positions 
at which the terminal velocity is reached for a given I. 
Although the distance between those terminal-velocity 
points and the tangent points is typically small, the non- 
circular motions and spiral shocks generate kpc-scale de- 
viations and discontinuities in the terminal-velocity lo- 
cus. Since those deviations happen at the positions of the 
spiral arms, they will generate larger errors at the vicin- 
ity of the arms, and will strongly affect the observer's 
view of the spiral structure of the model galaxy. 

For the outer galaxy (r > Rq), the usual procedure to 
determine the rotation curve involves looking for sources 
with independently k nown distance and m easuring their 
line-of-sight velocity ijBrand fc Blitd H993. for example). 
This procedure was simulated by assuming that the ob- 
server finds such a source at each point of the numerical 
grid outside the solar circle. It is assumed that the dis- 
tances to such sources are less reliable the farther they 
are from the observer. So, the circular velocity for the 
outer galaxy was taken to be, 

»^=i(x:2t+"°). w 

where the weights decrease linearly with the distance 
to the observer, and the summation is performed, at a 
given radius, over the azimuthal points excluding those 
within 7° of the galactic center and anti-center directions. 
Figure Ufa) shows the so obtained rotation curves, to- 
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Fig. 4. — a) The rotation curve given by the background potential (solid line), compared with the measured rotation curve. For the 
inner galaxy, both the rotation measured for positive (dotted line) and negative longitudes (dashed line) arc presented. The rotation curve 
corresponding to the outer galaxy is also shown (dashed-dotted line), b) The rotation curve given by the background potential (solid line), 
compared with the average of the north and south rotation curves (dotted line) and the mean azimuthal velocity of the gas in the simulation 
(dashed line). Notice that both the mean velocity and the mean rotation curve are above the background rotation curve for most of the 
radial domain. 



gether with the rotation consistent with the background 
gravitational potential. The northern rotation curve is 
lower than the southern rotation at ~ 3.5 kpc < r <~ 
5.5 kpc, while the opposite is true up to r = Rq. This 
behavior is sim i lar to the rotation curves reported by 
Blitz fc Spergel 1)1991^ . when scaled for R Q = 8 kpc. 

In order to try to recuperate the true (background) 
rotation, that should more closely trace the large scale 
mass distribution, the average both northern and south- 
ern rotation curves was taken. The result is compared 
with the background rotation in Figure Elb). Although 
the result is smoother and closer to the rotation consis- 
tent with the background potential, it is still systemat- 
ically higher (in agreement with the results reported by 
lSinhalll97Sfr . Another approach is to take the full vcloc- 
ity field and average the azimuthal velocity of the gas 
(Brand & Blitzl llQQST) . The result, also shown in Figure 
2jb) , is much closer to the background rotation, but it is 
still systematically larger. 

4. ERRORS IN THE KINEMATIC DISTANCE 

After adopting a rotation curve, and assuming that the 
gas follows circular orbits, the errors in the measured 
kinematic distances can be estimated by comparing the 
measured with the real distance in the simulation. 

In order to resolve the distance ambiguity for the in- 
ner galaxy, the usual procedure is to bracket the distance 
close enough as to place the object of interest on either 
side of the tangent point by l ooking at the gal actic latti- 
tude extension of the source ((Fish et alJl2f)08l). or using 
observed intermediate absorption features l|Watson et alJ 
120031 iSewilo et alJ 12004]) . For this investigation, I de- 
cided to cheat: I looked up on which side of the tangent 
point the gas parcel fell on, and chose the measured dis- 
tance accordingly. 

Figure |5fa) shows the error in measured distance with 



respect to the real distance in the model. Recalling Fig- 
ure the observer would determine different rotation 
curves for the northern and southern sides of the galac- 
tic center. Accordingly, in determining the kinematic dis- 
tance for Figure [SJa), the rotation curve used was that 
of the corresponding side of the galactic center. It is 
noticeable that although the errors are of the order of 
0.5 kpc in most of the galactic disk, they are significantly 
larger at the position s of the spiral arms (as hinted by 
iGomez fc Coxll2~004b|) . This fact has a special impact in 
studies of the spiral structure of the Galaxy that rely 
on kinematic distances, since it distorts the image the 
observer would generate (see M4.1J) . 

There is another significant feature in Figure EJ Al- 
though the terminal velocity does not really arises from 
the tangent point, the circular orbits assumption assigns 
gas observed near terminal velocity to that point. This 
fact generates a feature in the errors that corresponds 
to the locus of the tangent points. Again, the error is 
significant at the position of the spiral arms and would 
generate large errors in the determination of distances to 
objects that trace the spiral structure. 

The assumptions of circular orbits and different rota- 
tion curves for positive and negative longitudes are, of 
course, inconsistent. One solution is to fit a single ro- 
tation curve to both sides of the Galaxy. In order to 
test this method, the average of both rotation curves 
was taken and the equivalent of Figure Ela) was calcu- 
lated. The result was that the magnitude of the error in 
the kinematic distances was approximately the same, but 
the area with error > 0.5 kpc spanned a larger fraction 
of the disk. 

Suppose now that the imaginary observer somehow 
manages to obtain the large scale distribution of stellar 
mass in the model galaxy. This would allow the deriva- 
tion of the real rotation curve from the background ax- 
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Fig. 5. — Error in the measured kinematic distance (Ad) obtained under the assumption of circular orbits following (a) the measured 
rotation curve, and (b) the rotation given by the background potential. Sections of 7° around the galactic longitudes I = 0° and 180° were 
excluded. Although the error in most of the galactic disk is of the order of 0.5 kpc, it is significantly larger at the positions of the spiral 
arms. The sharp edges at the position of the tangent points are a consequence of the fact that the terminal velocities do not happen at 
those points. The errors in measured kinematic distances are larger when the real (background) rotation curve is used. 



isymmetric potential. If the observer now uses that real 
rotation to determine kinematic distances, even larger 
distance errors would be obtained, specially for the inner 
galaxy, as shown in Figure . This counter-intuitive 
fact arises because, at this point of the simulation, the 
gas has already adopted orbits that are not only influ- 
enced by the background potential, but also the large 
scale magnetic field (likely different from the field in the 
initial conditions) and the torques and resonances gener- 
ated by the spiral perturbation. Although the real rota- 
tion curve is consistent with the most important deter- 
minant of the gas rotation velocity (the background mass 
distribution), it does not include other influences in that 
velocity, while the "wrong" rotation curve determined 
from gaseous terminal velocities more closely reflects the 
real motion of the gas (recall Fig. 0] where the measured 
rotation curve is systematically above the true rotation) . 

Although intrinsically inconsistent, the two different 
measured rotation curves were used in the remaining of 
this investigation since that procedure leads to smaller 
distance errors. The results presented in the following 
section are even more notorious if the average or the real 
rotation curves are used. 

4.1. The Galaxy Distorted. 

Consider now that the imaginary observer tries to 
study the spiral structure of the galaxy he/she lives in. 
The procedure would consist of translating the longitude- 
velocity data obtained from a diffuse gas survey, for ex- 
ample, into a spatial distribution using the kinematic dis- 
tances that result from the assumption of circular orbits 
that follow the measured rotation curve. 1 The resulting 

1 In order to diminish spurious interpolation effects, each gas 
parcel was spread using a 2D gaussian weight function into a 3 X 3 
grid-cell region around the position corresponding to that parcel's 
galactic longitude and measured distance. 



map is shown in Figure Notice that the features de- 
scribed for Figure (namely the four spiral arms, the 
4 kpc high density ring, and the corotation low density 
ring) all but disappear, while new fictitious features, like 
the structure in the outer Galaxy, are formed as a con- 
sequence of the oscillations in the outer rotation curve. 
Also significant in this Figure are the regions where lit- 
tle or no gas is assigned by the mapping, namely the 
bands near the corotation circle and the quasi-triangular 
regions near the tangent point locus. (This nearly empty 
regions are significantly larger when the background or 
the mean rotation curves were used to determine the dis- 
tance to the observed gas parcel.) 

The imaginary observer would likely conclude that 
his/her home galaxy has 2 ill-defined spiral arms. If a 
logarithmic spiral model were forced, an ~ 11° pitch an- 
gle and a density contrast much stronger than that in 
the numerical model would be found. 

Another possibility for determining the distance to a 
gas parcel consists in comparing the line-of-sight velocity 
of the parcel with the predicted velocity obtained from 
some model for the galactic structure. For the numer- 
ical model described in given a galactic longitude, 
Figure is searched for the required velocity, and the 
corresponding distance is read out. 2 Although the same 
procedure to solve the ambiguity with respect to the tan- 
gent point was used, the non-circular motions introduce 
new distance ambiguities for certain longitude-velocity 
values (up to 11, although 3 is a more typical number). 
When these ambiguities appear, they happen close to 

2 A simple C language program that provides a distance given a 
galactic longitude and line-of-sight velocity value and uncertainty 
is available at http://www.astrosmo.unain.mx/~g.gomez/publica/ 
In that program, the resulting distance is given as a range, instead 
of a central value and uncertainty, since the velocity-distance map- 
ping makes the distance probability distribution neither uniform 
in the range, nor peaked around a central value. 
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each other, making their resolution difficult. So, when 
reconstructing the map of the galaxy, the gas density 
was equally split among these positions. 

The result is shown in Figured The new distance am- 
biguities still introduce spurious structure, like the split- 
ting of the spiral arms. Nevertheless, the number and 
position of the arms, the structure around the corota- 
tion radius and the lack of features in the outer galaxy 
are recovered. The imaginary observer would likely con- 
clude that his/her home galaxy has 4 arms with 9° 
and 12.5° pitch angles, although he/she would also find 
non-existing bridges and spurs. On the other hand, it 
should be considered that the imaginary observer does 
not see thermal nor turbulent line broadening. When 
these are considered, some of the new ambiguities will 
be swallowed into a distance range, effectively dissapear- 
ing. Therefore, some of the spurious structures will blend 
with real structures. So, the observer might get an im- 
age of the model galaxy closer to reality than Figure 
suggests. 

5. SUMMARY AND DISCUSSION 

The effect of the circular orbits assumption on our idea 
of the large scale structure of the Galaxy was explored. 
Since these errors might be quite large at the position 
of the spiral arms, the study of the spiral structure of 
the Galaxy and objects associated with it is particularly 
affected. By simulating the way an imaginary observer 
inside the model galaxy might try to infer the structure 
of the gaseous disk, it was found that the circular or- 
bits assumption destroys the spiral structure and creates 
spurious features in the measured distribution. 

The method of kinematic distances is a powerful 
one since it allows measurement of distances to diffuse 
sources and it is easily applicable to a large fraction of 
the galactic disk. Even if the measured rotation curve in- 
cludes deviations that do not reflect the true large scale 
mass distribution, Figure a) shows that the errors in 
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Fig. 6. — Re-mapping of the gas distribution resulting from the 
kinematic distances using the measured rotation curves in Figure 
[IJa) and assuming circular gas orbits. Notice the regions near the 
tangent points and the corotation circle, where little or no gas is 
mapped to. 
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Fig. 7. — Same as Figure © but using the full velocity field 
to recover the density distribution. Most of the characteristics of 
Figure are recovered, although some spurious structure appears 
due to the new distance ambiguities introduced by the non-circular 
motions. 



the distance are, in fact not very large for most of the 
galactic disk; in fact, the distance errors that arise from 
using the true rotation curve are larger. In both cases, 
however, the errors are quite large at the positions of 
the spiral arms. If we want to use this distance method 
for objects associated with the spiral structure, we need 
to conside r non-circular motions (as ha s been succesfully 
shown by iFoster fc MacWilliamil200fl for a set of H II 
regions and SNR). 

One possibility to achieve this is to try to determine 
the full velocity field of the galactic disk. But direct mea- 
sure of distances to the diffuse gas component is quite 
difficult (therefore the strength of the kinematic distance 
method). So, we need to use discreet objects and assume 
that they share their velocity with the diffuse component 
(Brand & Blitz 1993, for example; see also discussion in 
iMinn fc Greenbertall973|) . Yet another difficulty arises 
when tangential velocities and distances are required be- 
yond the solar neighborhood. 

Another approach at de termining the full veloci t y field 
is to model it. Recently, IFoster fc Mac Williams! (2006) 
used an analytical model of the density and velocity fields 
of the diffuse gas, with parameters for the model fitted 
to H I observations of the outer Galaxy. Despite the 
fact that their density and velocity models are not con- 
sistent in the hydrodynamics sense, and that the model 
do not include the dynamical effects of magnetic fields, 
they were able to add features of the Galaxy that are cur- 
rently difficult to incorporate to numerical models, like 
the disk's warp or the rolling motions associated with the 
spiral arms. Further numerical studies should allow the 
development of a more realistic analytical model. 

Instead of an analytic model, a numerical model was 
used in the present work to obtain density and velocity 
fields. Since the focus is in large scale velocity struc- 
tures, an eulerian code provides a good approach. Also, 
since the galactic magnetic field has been proved to 
be an important component of the total ISM pressure 
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ijBoulares fc Coxlfl99(]() . its effect in the gas dynamics is 
iikely to be important; therefore, a full MHD simulation 
was required. The large scale forcing is also trascen- 
dant; since the azimuthal shape of the spiral perturba- 
tion appears to have an influence on the gaseous response 
l)Franco et al.l l20?)3) . the usual sinusoidal perturbation 
was deemed too simplistic and a self-consistent model for 
the perturbing arms was chosen. At the present time, 
the galactic warp an d the vertical motions associated 
with the spiral arms ijGomez fc Coxll2004albl) could not 
be considered at the necessary resolution. 

In this work, it has been shown that it is possible to 
recover most of the gaseous structure of the galactic disk 
using kinematic distances, as long as the full velocity 
field is considered. Nevertheless, applying these results 
to the Milky Way is a whole new issue, since obtaining 
the full velocity field is not trivial. For the procedure 



used here, how close the numerical simulation is to the 
real Galaxy remains the weak point of this approach. 
The computation cost of a realistic enough simulation 
is still too high to allow a parameter fitting analysis. 
So, the remaining question is if the velocity field that 
results from the simulation yields a determination of the 
distance to a given object, or only an estimation of the 
distance error. The answer to that question is left to the 
reader's criterion. 
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